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Abstract 

A study of collisionless external shocks in gamma-ray bursts is presented. The 
shock structure, electromagnetic fields, and process of electron acceleration are assessed 
by performing a self-consistent 3D particle-in-cell (PIC) simulation. In accordance with 
hydrodynamic shock systems, the shock consists of a reverse shock (RS) and forward 
shock (FS) separated by a contact discontinuity (CD). The development and structure 
are controlled by the ion Weibel instability. The ion filaments are sources of strong 
transverse electromagnetic fields at both sides of the double shock structure over a 
length of 30 - 100 ion skin depths. Electrons are heated up to a maximum energy 
tele ~ where e is the energy normalized to the total incoming energy. Jet electrons 
are trapped in the RS transition region due to the presence of an ambipolar electric field 
and reflection by the strong transverse magnetic fields in the shocked region. In a process 
similar to shock surfing acceleration (SSA) for ions, electrons experience drift motion 
and acceleration by ion filament transverse electric fields in the plane perpendicular 
to the shock propagation direction. Ultimately accelerated jet electrons are convected 
back into the upstream. 

Subject headings: acceleration of particles — galaxies: jets — gamma rays: bursts — 
magnetic fields — plasmas — shock waves 
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Introduction 


Relativistic jets are present in a variety of high-energy astrophysical systems, e.g., pulsar 
wind nebulae (PWNe), gamma-ray bursts (GRBs), and active galactic nuclei (AGNs). The jets 


(Bromberff et a 

20111 

(Ramirez et al. 

2002 

), 


powerful internal jets and supernovae explosions. Non-thermal emission associated with the rel¬ 
ativistic jets is believed to be generated by the accelerated electrons via synchrotron or inverse 
Compton mechanisms. Diffusive shock ac celeration (PSA) is the m ost applicable theory for the 
acceleration of particles in shock systems (|Blandford &: Eichleiill987l b In this process, particle ac¬ 
celeration takes place as the result of particle trapping and oscillating at the shock front. Even 
though DSA is successful for elucidating ion acceleration, acceleration of the electrons cannot be di¬ 
rectly related to this mechanism. To experience DSA a thermal electron requires crossing the shock 
front multiple times. However, a thermal electron is closely tied to the magnetic field lines because 
of its small Larmor radius. Thus with no strong pre-acceleration process, a thermal electron is 
convected downstream without feeling any serious DSA. The lack of DSA for t hermal electrons is 


Guo et al.ll2ni4l ). 


often referred to as th e electron injection problem (jHoshino fc Shimada 
2009 


2002; 


Amano fc Hoshino 


PIG simulations can shed light on the mechanism of electron acceleration in collisionless shocks. 
In the high Mach number (the ratio of the upstream bulk flow speed to the sound speed in the 


energy electrons into the DSA ( 

Dieckmann et al. 

2000; 

Hoshino 

2001; 

McClements et a 

2001 

Hoshino & Shimada 

2002: 

Schmitz et al. 

2002: 

Amano & Hoshino 

2009; 

Matsumoto et al. 

20 d) 


In the SSA process predominantly electrostatic flu ctuations are excited at the leading edge of the 
shock via the Buneman instability (lBunemanlll958l ). due to the interaction between incoming elec¬ 
trons and reflected ions. The incoming electrons are then trapped in the electrostatic potential and 
efficiently accelerated by the shock convective electric field (i.e., the electric field resulting from the 
motion of the magnetized upstream region toward the shock, Eq = —/3o x Bo). 

To date, electron acceleration in unmagnetized (Bq = 0) electron-ion shocks has been poorly 
investigated. The injection process is anticipated to be different from magnetized shocks because 
there is no convective electric field to accelerate the electrons. In the present work, a 3D PIG 
simulation is used to study the shock structure, and the process of electron acceleration associated 
with a relativistic jet propagating into an unmagnetized plasma. This setup is the most applicable 
model for external GRBs shocks. The spati al and temporal scales of the present simulation are much 


large r than those used in previous studies (jNishikawa et al 


2003 


Ardaneh et al.ll2014l : iGhoi et al 


2014l i in order to allow us to observe efficient particle acceleration. The present paper aims to 
address the following questions related to an unmagnetized double shock system: (i) What are the 
field structures responsible for the processes of electron heating, and acceleration?, (ii) Where do 
these processes mainly take place?, (hi) What is the resulting associated electron energy spectrum?. 
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and (iv) What are the principal mechanisms responsible for electron heating, and acceleration? 

In this paper the simulation method and parameters setup are summarized in § 2. The main 
simulation results are discussed in § 3. Finally, we conclude in § 4. 


2. Simulation method and parameter setup 


In recent works, in order to generate a shock, a relativistic plasma stream is injected from 
one end of the computatio nal grid and reflecte d from a rigid wall at the opposite end (for in- 
stance, a ID simul a .tion by Hoshino k, Shimada (2002), 2D simulations by Spitkovskv ( 2008al f^: 
Amano &: Hoshino (2009); Martins et al. ( 20091 '). and a 3D simulation bv ICuo et ahl (j2014 ')l. This 
method resembles the collision of two identical counter-streaming beams and reduces by one-half 
the number of calculations. This approach simulates only a moving FS. In the present work we 
inject a particle jet into an ambient plasma, so that a double shock structure is fully captured. 
With this setup the jet-to-ambient density ratio can be changed and by assuming a density ratio 
greater than one, the shock formation process can be properly handled by smaller scale (temporal 
and spatial) simulations. Deceleration of the jet flow by the ambient plasma results in a CD (the 
CD is the location where the electromagnetic field and the velocity of jet and ambient plasmas are 
similar and the density changes). Two shocks propagate away from the CD into the jet and ambient 
upstreams (in the CD frame). Hence, two shocks and one CD split up the jet and ambient plasma 
into four regions: (1) unshocked ambient, (2) shocked ambient, (3) shocked jet, and (4) unshocked 
jet. Hereafter, subscripts 1, 2, 3, and 4 refer to the unshocked ambient, shocked ambient, shocked 
jet, and unshocked jet, respectively. Quantities with a single index Q\ denote the value of quantities 
Q in region i in frame i and quantities with double indices Qy denote the value of quantities Q in 
region i as seen in frame j. 


The present simulatio n is performed using an adopted version of the TRISTAN code (jBuneman 


1993 


Niemiec et al.l 120081 1. a massively par allel, fully relativistic , PIC code that can be used for 


many application s in astrophysical plasmas (jNishikawa et al.ll2003l . l2005l . [2009j; lArdaneh et al.ll2014l : 
Choi et al.l 120141 ). The simulation is carried out in the ambient frame with a box including 8005 
nodes in the x-direction, and 245 nodes in the y and 2 ;-directions. There are 6 particles per cell per 
species for the ambient plasma (~ 3 billion particles per species) and the jet-to-ambient density 
ratio used is 1.7. The ion-to-electron mass ratio used is 16. The system is numerically resolved 
with 5Ax per electron skin depth in each direction and AtUpe = 0.01, where Ax and At are the 
grid size and the time step, respectively. The thermal speed of ambient electrons and ions is 0.05c 
and 0.01c, respectively. In the simulation, the jet fills the whole computational box in the y 2 ;-plane 
and is injected at x = 25Ax in the positive x-direction. The injected jet bulk speed is = 0.9950, 
741 = (1 — = 10, the Mach number M = /?4i//3si = 60, and the jet electrons and ions 

have a thermal speed of 0.014c and 0.003c, respectively. The surfaces at x = Xmin and Xmax are 
rigid reflecting boundaries for the ambient particles, while they are open boundaries for the jet 
particles. These surfaces are radiating boundaries for the fields based on Lindman’s method as 
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explained in iBunemanI (jl993l i. Periodic boundary conditions are applied at all other boundaries 
for both particles and fields. In the following, time is normalized to space to the electron skin 
depth, Ace = c/cjpe, particle momentum for species i to the corresponding mic (e: electron and i: 
ion), and density to the upstream ambient density, ni. The dimensionless form of a variable, such 
as X, is shown by x*. The energy is normalized t o the total incoming energ y density, ^ 
and is denoted by e, according to the notation of iMedvedev fc LoebI (jl999l ). 


3. Simulation Results 

The present simulation includes growth of the Weibel instability, and the generation of elec¬ 
tromagnetic fields which decelerate the jet flow and consequently form a double shock structure. 
At late times the particles are effectively heated, and accelerated. In the following this scenario is 
clarified in more detail. 


3.1. Shock structure 


When the jet with density ratio of n 4 i/ni = 1.7 and 741 = 10 interacts with the ambi¬ 
ent plasma, the total density initially increases to 2.7ni near the jet front. In this interac¬ 
tion, the distribution of particles is extremely anisotropic and this distri bution is unstab le to 
several plasma instabilities, such as the electrostatic Buneman instability (IBunemanI Il958l ) and 
the electromagnetic Weibel instability ([Weibell Il959l j. In relativistic shocks, the Weibel insta¬ 


bility has the large s t growth rate and will dom inate the interaction (|Medvedev &: LoebI Il999l : 
Califano et al.l 12002: iHededal &: Nishikawal l2005l b Time evolution of the Weibel instability was 
analyzed in lArdaneh et al.l (|2014| j and we summarize the results here. Jet electrons are deflected 
when interacting with ambient particles to form electron current filaments. Subsequently magnetic 
fields are amplified principally as a result of mutual attraction between the electron filaments and 
their merging to larger scale filaments. When the magnetic fields become strong enough to deflect 
the heavier ions, the ions begin to participate in the instability. The ion current filaments then grow 
with a smaller growth rate related to the electron-ion mass ratio and to transverse heating. As time 
evolves, the Weibel instability induced electromagnetic fields facilitate momentum transfer between 
the jet and ambient plasma. A double shock structure forms when the relative velocity between 
the jet and ambient plasmas exceeds the sound speed in the ambient plasma and the magneto-sonic 
speed i n the jet, and the pressur e in the shocked region exceeds the pressure in the unshocked 
region (jZhang &: Kobavashil l2005l b A first density compression appears due to the deflection of 
jet electrons and ambient ions which we define as the RS. This density compression grows due to 
amplification of the Weibel instability at the expense of the streaming ions. Time stacked plots 
of the transversely averaged (in the yz-plane) total ion density as a function of axial distance are 
shown in Fig. [T^. The time range is t* = 20 — 540 with an interval of At* = 20. One can readily 
identify a RS propagating in the positive x-direction with /3rsi ~ 0.70. The shape and peak value 
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of the total density corresponding to the RS are almost constant in time {n-n/nn « 2.35, Fig. [TJ)), 
and the RS is fully developed. On the other hand, a density compression appears mainly in the 
ambient plasma at late times {t* ~ 300) that we identify as the FS. The compression ratio rises 
with time until reaches about n 2 i/ni ~ 6 at the end of the simulation (Fig. Hb)- The FS structure 
moves at a speed /3 fsi ~ 0.89 in the positive x-direction. The CD moves in the positive x-direction 
at a speed /3cdi ~ 0.80. The shocked region (340 < x* < 430 at t* = 500) between the RS and FS 
separates the jet and ambient upstreams (Fig. [Tb). 




Fig. 1.— Panel (a) shows stacked profiles of the transversely averaged ion density (total 
jet+ambient) from t* = 20 up to 540 with an interval of At* = 20 and panel (b) represents 
the ion density profile at t* = 500. Solid, dashed-dot, and dotted black lines in panel (a) indicate 
the RS, FS, and jet front, respectively. In panel (b) the shocked region is shown as orange, yel¬ 
low represents the shock transition region (left for the RS and right for the FS), white shows the 
unshocked regions (jet or ambient), and the dashed blue line shows the CD. 


The results of present simulation are quantitatively validated by the theoretical calculations 
described in the Appendix. The speed of the CD in the ambient frame from Eq. (IA7p i is /3 cdi = 0.8. 
Combining Eq. (|A3h l and (|A3b l , the ES speed can be determined. Eor a shocked ambient adiabatic 
index between 4/3 < F < 5/3 and 7cdi = 1-67 (corresponding to /IcDi = 0.8) the FS speed will 
lie between 0.85 < /Ipsi < 0-90 (1.91 < ypsi < 2.23) where upper and lower limits correspond to 
the upper and lower limits of F, respectively. From Eq. (IA8I) . the jump condition at the ES lies 
between 5.66ni < n 2 < 9.66ni, where n 2 is the shocked ambient density in the CD frame and the 
lower and upper limits correspond to the upper and lower limits of f, respectively. Measured in the 
ambient frame the density compression would lie between 9.4ni < n 2 i < IG.Oni. In the jet frame 
the CD moves with speed /3 cd 4 = —/34CD = —0.96 (Eq. (IA5h D and the corresponding Lorentz 
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factor is 7 cd 4 = 3.41. By 1 —>■ 4, 2 —>• 3 and 734 ^ 7CD4 transformations in Eqs. ()A3p . we can 
calculate the speed and the corresponding Lorentz factor of the RS in the jet frame. In this frame 
we have 0.97 < —/ 3 r,s 4 < 0-98 and 4.35 < 7 rs 4 < 5.77, where upper and lower limits correspond to 
upper r = 5/3 and lower F = 4/3 limits for the shocked jet adiabatic index, respectively. The RS 
speed in the ambient frame (Eq. (lAlOM lies between 0.51 < /3 rsi < 0.70. The density jump at the 
RS region can be obtained from Eq. ()A9p where 744 = 7441 / 741 . Hence the density compression of 
shocked jet is between I.O 7441 < 743 < 1 . 667441 according to the upper and lower limits of T. In the 
ambient frame the shocked jet density lies between 1 . 77441 < 7431 < 2.7?44i. Although the density 
jump for the shocked ambient is about a factor of ~ 1.5 smaller than theoretically predicted for a 
fully developed ES, the simulation results for the RS, CD, and ES speed and also for the density 
jump at the RS are in the ranges provided by the theoretical analysis for a fully developed shock 
system. 


Shown in Eig. [2^ and [2]3 are the z-component of the electric field, and y-component 
of the magnetic field, By, respectively. Panels [2b-l2^ show the transversally averaged (in the yz- 
plane) electric and magnetic field components, [Ex(Rx), T^y(Rz),T'z(Ry)]- The energy distribution 
(total of jet+ambient) and average energy along the x-direction for electron and ion species are 
shown in Eigs. [2fj2l. All panels are at t* = 500. Where high-speed jet particles interact with 
the ambient medium (behind the RS at x* < 340) or scattered ambient particles blend with 
the upstream ambient (in front of the ES at 430 < x*), particles distribution becomes strongly 
anisotropic. Anisotropies result in the Weibel instability which generates current filaments in these 
regions with currents in the x-direction. According to Ampere’s law, these current filaments are 
encircled by transverse magnetic fields, and we see that (Rx) = 0 in Eig. [2j:. The transverse electric 
fields are related to the magnetic fields via /3e:i x E = B where /3e:i is the velocity of the electron 
(ion) carrier. The carriers move roughly at the speed of light in the x-direction, /3e:i — /3e:ix — 1- 
Therefore, the transverse electric field components are Ey = B^, and E^ = —By, as are observed 
in the simulation results for [Ey{Bz), Ez{By)] (Eigs. [2^-[2^). Additionally, there is a longitudinal 
ambipolar electric field within the RS transition region, 140 < x* < 340 for t* = 500 (Eig. ^). 
This electric field is generated by the density gradient and different mobilities of electron s and 
ions ( Enrslimd fc Shnnk 1970 : Enrslimd &: Ereidberg 1971 : Hoshino 2001 : Choi et al. 2014). The 
magnetic fields act to isotropize the momentum distribution, while the electric fields function to 
thermalize, and accelerate the particles afterwards. In Eigs. [2]F-|2}, the shocked region lies between 
X* = 340 and x* = 430. Within the RS transition region (140 < x* < 340) jet electrons are 
trapped by the ambipolar electric field and effectively accelerated up to 7 e = 200 by the transverse 
electric helds (Eigs. [2f and [2^). A tenuous population of these electrons convect upstream due to 
reflection by the magnetic fields in the shocked region (ellipse in Eig. [2f). On the other hand, jet 
ions are slowed in the RS transition region (140 < x* < 340) by 40% from the initial Lorentz factor 
7 i = 10, due to the effect of the ambipolar electric field. In the shocked region, jet electrons have 
been fully thermalized and are well merged with the thermalized ambient electrons. Thus, only 
a single electron population is present in the hot shocked region (Eig. [2f). On the other hand, 
the kinetic energy of jet ions is transferred to the heating of ambient particles by means of the 
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electromagnetic fields generated by the ion Weibel instability (Fig. [2)r). Full thermalization of the 
two ion populations (jet and ambient) has not yet occured (demands a longer simulation time), i.e., 
the two populations are distinguishable in Fig. [2)i. Electrons located in the FS transition region 
(430 < X* < 500), predominantly ambient electrons, also undergo the Weibel instability and are 
thermalized by the jet upstream kinetic energy. In this region, penetrating jet ions interact with 
ambient particles and are slowed down gradually by 50% from the jet front Lorentz factor 71 = 10 
to a minimum value of 71 = 5 (Figs. [2)i and[2j). 



Fig. 2.— Structure of the relativistic double shock at t* = 500. Panels (a) and (b) show the z- 
component of the electric field and the y-component of the magnetic field By at y* = 24. Panels 
(c)-(e) show the transversally averaged (in the yz-plane) field components. Total of the jet (blue) 
+ ambient (red) particle energy distribution, and average energy along x-direction are presented 
in panels (f) and (g) for electrons, and (h) and (i) for ions. The shocked region is identified by 
orange coloring or between vertical orange lines (as shown in Fig. [1]). The ellipse in panel (f) 
shows the high-energy electrons reflected into the upstream. In (f) and (h), 1.2 x 10® particles 
are randomly selected. Due to the very large number of particles in the simulation, only a part 
of them is represented. Particles of each population are selected randomly so that the respective 
distribution function is not affected. 


Ion current filaments generated by the ion Weibel instability play a crucial role in establishing 
the shock transition, electron heating, and acceleration in an unmagnetized electron-ion plasma 

































































( Frederiksen et al.ll2004 : Hededal et al. 2004 : Spitkovskv 2008a). In Figs. [3l the longitudinal current 
density, Jx = Jix + Jex, through a transverse cross section is shown at x* = 320, in the RS 
transition region (Fig. [3^ and [3b) and at x* = 480, in the FS transition region (Fig. [3b> and[3]l) for 
t* = 500. Positive (red-white colors) and negative (green-blue colors) represent the ion and electron 
contribution to the total current, respectively. As one can see, longitudinal current filaments are 
surrounded by approximately azimuthal magnetic fields (illustrated by the arrows in Fig. [3b and 
[Sb). Electric fields are perpendicular to the magnetic fields and the associated arrows point in the 
radial direction (Fig. [3b and[3bl). In the RS transition region, magnetic fields are predominantly 
due to the ion filaments (counter-clockwise arrows) and electrons act to Debye shield these filaments 
([Frederiksen et al.[[2004b . However, there are some electron filaments (clockwise arrows) within the 
FS transition region. Furthermore, electromagnetic fields in the RS transition region are stronger 
than those in the FS transition region. The transverse size of the filamentary structures is on the 
order of several electron Larmor radii (R* ~ 3 — 5). Perpendicular electric and magnetic fields lead 
to the E X B drift motion of electrons parallel to the shock propagation direction (x-direction). 
During this motion, the electrons are effectively heated. 


3.2. Acceleration mechanism 


One of the key factors in generating non-thermal radiation from GRB afterglows is a non- 
thermal, high-energy electron population. This population can be seen in the electron spectrum, 
where a pure 3D Maxwell-Jiittner distribution does not account for the high energies. In fact, a 
more complex spectrum is anticipated as a result of electron acceleration. Shown in Fig. [4|is the 
electron spectrum from time t* = 420 up to 500 with an interval of At* = 20 taken at 200Ace 
in the RS transition region where most of the accelerated electrons are located. The electron 
spectrum consists of a main thermal peak around p* ~ 14, and a high-energy tail which extends 
to p* PS 400 at t* = 500. The dashed line shows a 3D Maxwell-Jiittner distribution with f{p*) = 
p*"^/T*K2{1/T*) exp (—Y'T+"p^/r*) where K 2 is the modified Bessel function of the second kind, 
and T* = K-^T/nief?- We set T* = 5, in accordance with T* = (741 — l)n 4 i/ 3 ni = 5 which is 
expected if the upstream jet kinetic energy is completely converted into internal kinetic energy of 
the jet and ambient particles. The dashed-dotted line shows a power-law fit to the non-thermal, 
high-energy electron population. The power-law begins around p* = 100 and extends to p* = 360 


Kirk et al. 

2000 : 

Achterbere et al. 

2001 ) 


The mechanism responsible for the non-thermal tail is analyzed by tracing the trajectories 
of some of the electrons located in the RS transition region at t* = 420. A t first, the electrons 
are heated with heating principally associated with the ion current filaments ([Hededal et al.[[20041 : 
Medvedev[ l2006l: [Soitkovsk-vi [2008a[ i . As discussed previously, when a relativistic jet interacts with 
the ambient medium, current hlaments are formed by the Weibel instability. The ion filaments 
are sources of strong electromagnetic fields in the transition regions of the RS and FS, while the 
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Fig. 3.— Longitudinal current density through a transverse cross section at x* = 320 (hrst column) 
and X* = 480 (second column) for t* = 500. The arrows represent the transverse magnetic (first 
row), and electric (second row) fields. 


electrons provide Debye shielding of the ion filaments. The electromagnetic fields in the vicinity of 
the ion current filaments can be seen in Figs. [3j The magnetic fields are predominantly azimuthal 
while the electrostatic fields are in the radial direction. The presence of perpendicular electric 
and magnetic fields results in a drift motion of the electrons with Ud = |E x parallel to 

the shock propagation direction (x-direction). At distances less than a Debye length, an electron 
initially moves in a direction opposite to the electric field. Due to the velocity v± thus acquired, the 
magnetic part of the Lorentz force produces a motion perpendicular to the electric and magnetic 
fields which instantly bends the electron trajectory. During the first-half gyration, an electron is 
accelerated at the expense of the potential energy stored in the transversal electric field (A-B in Fig. 
[5^ and Fig. [ 5 ) 3 ). On the other hand, an electron in its second-half gyration does work on the electric 
field and its kinetic energy will be transferred to the electric field energy (B-C in Fig. [5^ and Fig. 
[ 5 ) 3 ). The growth of the ion Weibel instability increases the electric charge of the ion filaments with 
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Fig. 4. — The electron spectrum taken at 200Ace within the RS transition region from t* = 420 
(leftmost solid black line) up to 500 (rightmost solid red line) with an interval of At* = 20 (200Ace 
is calculated relative to the RS position at each time). The dashed line shows a 3D Maxwell- 
Jiittner distribution. Dashed-dotted line shows a power-law fit to the non-thermal component in 
the electron energy spectrum at the latest time. 


time, hence in the next gyration (C-D-E in Fig. [5^ and Fig. [Sb) an electron experiences deeper 
electric potentials that enhance the amplitude of energy oscillations. The maximum attainable 
energy for an electron during the heating stage can be estimated analytically. Approximately, an 
electron moving toward an ion current filament gains energy rigie ~ eEAr ~ eBAr (/3xi — !)• The 
maximum radial distance that an electron can travel is about half the distance between the filaments 
(measured from the filaments axes), Ar ~ c/wpi^r, where Wpi^r = (47rnie^/mi7i)^/^ is the relativistic 
ion skin depth. Hence, the electron energy density is Ceie ~ yTb (normalized to the total incoming 
energy). In the present simulation, the maximum of yTb for 140 < x* < 340, where electrons 
are found at high energies, is about 0.06 (calculated based on the magnetic field components in 
Figs. [Ib-Eb and normalized to the total incoming energy = 170). Therefore, the 

maximum Lorentz factor of electrons during the heating stage is about 102, in agreement with the 
3D Maxwell-Jiittner distribution dashed line shown in Fig. 01 


Acceleration of electrons within the RS transition region (Fig. [2]f) is associated with an ambipo- 
lar large-scale elect ric field induced in the negative x-direction due to the inertia difference between 
electrons and ions (jHoshino Il200lh . The basic idea behind the electron acceleration process is il¬ 
lustrated in Fig. [6b. In the figure the yellow region represents a distance in front of the ambipolar 
electric field and behind the RS (such as 250 < x* < 340 in Fig. [2]F) in which electrons can oscillate. 
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Fig. 5.— Panel (a) shows the 3D trajectory of an electron (selected at = 420) during heating 
over a time interval At* = 80 and panel (b) shows the electron energy over time. The color of the 
electron trajectory expresses its energy based on the color in panel (b). The longitudinal electron 
position is calculated relative to the RS front position (x* — /3rsi x t*). 


Consider an incoming electron that enters this region with a small /3ex- The electron moves in the 
positive x-direction (opposite the ambipolar electric field direction), and will be reflected by the 
strong magnetic fields {By,Bz) in the shocked region, and then this process may repeat. During 
this process, the electron experiences a drift motion in the yz-plane parallel to the RS plane (Figs. 

andEla). As a result, the electron is efficiently accelerated by the transversal electric fields (Fig. 
[6t) because for |E_l| = |B_l| (Figs. [2^-[2^) and 0 < /3ex < 1, leE^l > |e/3ex x B±|- After multiple 
reflections, the magnetic force experienced by the electron, e{j3eyB^ — PezBy) becomes larger than 
the electric force, eE^. At this point the electron will be convected in the negative x-direction and 
again be trapped, but now in the region of the ambipolar electric field (like 140 < x* < 250 in 
Fig. [2]F). The electron undergoes the same acceleration process in this region until it is convected 
into the upstream. The resulting small high-energy electron population is indicated by the ellipse 
in Fig. [2f. This acceleration pro cess resembles the SSA process for ions (|Sagdeev &: ShaDirdll974l : 


Zank et al.lll996l : iLee et al.lll996l i . This efficient accel eration pr o cess f o r electrons in the shock tran - 


sition region is indicated in I D PIC simulations bv iHoshino I (1200111: iHoshino fc Shimadal (j2002l i 


and in 2D PIC simulations by Amano &: Hoshino ( 2009l k 


Matsumoto et al 


(|20li)), 


where a series 

of large amplitude electrostatic solitary waves (ESWs) were induced by the Buneman instability 
resulting from the interaction between reflected ions and incoming electrons. Instead of the up¬ 
stream convective electric field in the SSA, in the present acceleration process there are transverse 
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electric fields generated by the Weibel instability that effectively accelerate the electrons. Although 
previous work has shown that at late times both electrons and ions are accelerated by the DSA 
(|Sironi et al.ll2013l i. in the present simulation, electrons are accelerated to the highest energies via 
the SSA process. 



Fig. 6.— Panel (a) represents the schematic diagram of the electron acceleration process, panel 
(b) is the 3D trajectory of an electron (selected at tg = 420) during acceleration in a time interval 
At* = 80, and in panel (c) the electron energy over time is illustrated. The colors of the electron 
trajectory reflect its energy according to the color in panel (c). The longitudinal electron position 
is calculated relative to the RS front position (x* — /3rsi x t*). 


4. Summary and conclusions 

In this work, we have studied the physics of electron acceleration in a relativistic collisionless 
Weibel instability mediated shock by means of a 3D PIC plasma simulation. In our simulation, an 
unmagnetized jet plasma with a 3D drifting Maxwell-Jiittner distribution and temperature Tg = 
is injected into an ambient unmagnetized medium. This process results in a double shock structure 















- 13 - 


which consists of three regions. 1) A pre-shock acceleration region (behind the RS), consisting 
of the incoming jet and ambient particles and with large particle distribution anisotropy. In this 
region, strong electromagnetic fields due to the Weibel instability, and an ambipolar electric field 
due to the density gradient and the different mobilities of electrons and ions are induced. 2) A 
fore-shock region (in front of the FS) where the magnetic fields are relatively weak but still strong 
enough to be considered for some GRB afterglow radiation. 3) A hot shocked region where the 
magnetic field energy density is a few percent of equipartition with the kinetic energy density. This 
region is clumpy, with large structures that are on the order of 30 ion skin depths in length. In 
the present simulation, the jet RS forms first. The ambient medium FS forms at later times. In 
this simulation, the RS is fully developed and the shock structure and jump conditions confirm 
theoretical predictions obtained from an analytical treatment. The FS still developing at the end 
of the simulation and has not yet reached the theoretically predicted jump conditions. 


Ion current filaments are the origin of deeply penetrating electromagnetic field structures. 
Heating of electrons is directly tied to the ion filaments. This links electron heating closely with 
the micro-structure of relativistic collisionless shocks. With no restriction on the stage of shock 
formation, the maximum electron energy density is Cgie ~ -^/ib during the electron heating stage. 
This is in excellent agreement wi th values of Cgie/y^ derived from afterglow data from 10 GRBs 


(jPanaitescu 


2001 


Medvedevll2006l f . Thus, there are two equipartition parameters in the GRB shock; 


the magnetic field energy density , and the energy density of the electrons, Ceie, normalized to 
the total energy density of the shock, respectively. In the current simulation, the electrons are 
heated by the Lorentz factor 7 e = 100. 


We have shown that non-thermal electrons with a power-law index of 2.2 are generated behind 
the RS. Concerning a relativistic jet propagating into an unmagnetized plasma, this non-thermal 
electron population is observed for the first time to the knowledge of the authors. Although the 
particles are mostly accelerated in the shock-dominated systems via DSA, tracing the electron 
trajectories reveals that the mechanism responsible for electron acceleration differs significantly 
from DSA. Acceleration of electrons is due to the presence of a strong ambipolar electric field 
behind the RS. The incoming jet electrons move opposite to the direction of the ambipolar electric 
field (x-direction in the present work). The electrons are then reflected back by the strong transverse 
magnetic fields in the shocked region. This results in an electron drift motion in a plane parallel 
to the RS plane while they are efficiently accelerated by the transverse electric fields. When 
the magnetic force becomes large enough to overcome the force of the ambipolar electric field, the 
electron are convected into the jet upstream of the RS. Therefore, a high-energy tenuous population 
of electron exists in the jet upstream. The acceleration mechanism is similar to the SSA for ions, 
where transverse electric fields generated by the Weibel instability play the role of the upstream 
convective electric field in the SSA. 


The authors thank the referees for insightful comments and suggestions that led to correct 
interpretation of the acceleration observed in the present simulation. The authors thank Philip 
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AST-0908040, NNX12AH06G, NNX13AP21G, and NNX13AP14G. The simulation presented here 
was performed using the supercomputer of AGGMS, Institute for Information Management and 
Communication, Kyoto University. 


A. Jump conditions for a relativistic 90°shock 


We review first the MHD analysis of a relativistic 90°shock propagating into an unmagne¬ 
tized plasma, and then extend the analysis to double shock structure. In the following, indices 
1, 2, and s refer to the upstream, downstream, and shock frame, respectively. The jump con- 
ditions for 90°shocks are solutions to the following one-dimensional MHD conservation equations 
( Kennel &: Coroniti 1984 1: 


[1 + 


[1 -I- cri(l - = 72 s /^2 

(1 - + ^ = + (—)^ 
Is Pis P2s 


2/3-2 


(Ala) 

(Alb) 


where cri = Bf ^/is the magnetization, with Bis as the transverse magnetic field, /i* as 
the dimensionless specific enthalpy, and Y as the ratio of shock frame magnetic fields or density 
( Kennel fc Goroniti 19841 ). where /i* and Y are defined by 


Y = — 

Bis 


* 111 1 

Pi = 1 + (—- 

ni'Ti-l 


P* 


N2s ^ l2sPls ^ 72sn2 
A^is 7 isP2s 7isni 


(A2a) 

(A2b) 


Here P* denotes the dimensionless thermal pressure, P* = P/nimcP, and f is the adiabatic index. 
In our simulations the upstream flow is considered a cold plasma, i.e., P^ = 0, so that = 1. 
Solving Eq. (jAlh l for and inserting t he resulting expression into Eq. (lAlb l leads to the following 
equations for an unmagnetized plasma (jZhang &: Kobavashill2005l l: 


*2 ^ (721 - i)(r 72 i +1)^ 
f(2-f)(72i-l) + 2 ’ 
2 ^ (721 + i)[r( 72 i -1) +1]^ 
f(2-f)(72i-l) + 2 


(A3a) 

(A3b) 


The corresponding jump condition can then be determined by Eq. ()A2b l. Now we consider an 
unmagnetized jet with Lorentz factor 741 being decelerated by an ambient medium with density ni. 
Eor such system, indices 1, 2, 3, and 4 refer to the unshocked ambient, shocked ambient, shocked 
jet, and unshocked jet, respectively. Relative to the CD, the above formalism can easily include 
the shocks associated with jet-ambient interactions. The speed of the CD can be determined by 
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applying Eq. (jAlb i at the interface to the ambient and jet plasmas (jNishikawa et al.ll2009l i. Hence, 
in the CD frame we have 


/^i7iCd/^1CD + A* — -/^474Cd/^4CD + ■ 

ni 

The jet speed measured in the CD frame and the associated Lorentz factor are given by 

^41 - ,dlCD 

- 1 - ftlAcD ’ 

74CD = (1 - ■ 


(A4) 

(A5a) 

(A5b) 


In our simulations, the jet and ambient medium are initially cold plasmas. Hence Eq. (IA4I) reduces 
to 

l^r7lCD/?lCD = 74 Cd/?4CD ) (^6) 

where is the jet-to-ambient enthalpy ratio. Making use of Eqs. ()A5p . after some algebra Eq. 
(jA 6 p gives 

1/2 

/5iCD = 

741/Xr^ + 1 

which is the speed of ambient medium towards the CD, and thus the speed of the CD through 
the ambient medium. In the shocked ambient, once /3icd (7icd) is determined from Eq. (|A7p . by 
721 —>■ 7iCD Eqs. (|A3p can be solved to find the speed of the shocked ambient (FS) in the ambient 
frame, /3fsi = —Pis- The jump condition for the FS is 


^ _ fyicD + 1 
ni f - 1 


In the shocked jet (RS), making use of Eq. (jA5h i allows the speed of CD in the jet frame to be 
calculated / 3 cd 4 = —PiCD- Under 1 —)• 4, 2 ^ 3 and 734 —>■ 74 CD transformations in Eqs. (|A3p . we 
can calculate the speed of the RS in the jet frame, / 3 rs 4 = —Pis- In this region the jump condition 
becomes 


^ _ r7cD4 + 1 
n 4 f - 1 


(A9) 


Finally, the speed of the RS in the ambient frame is given by 


Prsi 


Pii - Pnsi 
1 — PuPrsa 


(AlO) 
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